Open Mushrooms: Stickiness revisited 
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We investigate mushroom billiards, a class of dynamical systems with sharply divided phase space. 
For typical values of the control parameter of the system p, an infinite number of marginally unstable 
periodic orbits (MUPOs) exist making the system sticky in the sense that unstable orbits approach 
regular regions in phase space and thus exhibit regular behaviour for long periods of time. The 
problem of finding these MUPOs is expressed as the well known problem of finding optimal rational 
approximations of a real number, subject to some system-specific constraints. By introducing a 
generalized mushroom and using properties of continued fractions, we describe a zero measure set 
of control parameter values p G (0, 1) for which all MUPOs are destroyed and therefore the system 
is less sticky. The open mushroom (billiard with a hole) is then considered in order to quantify 
the stickiness exhibited and exact leading order expressions for the algebraic decay of the survival 
probability function Pit) are calculated for mushrooms with triangular and rectangular stems. 

I. INTRODUCTION 

Billiards [JJ are systems in which a particle alternates between motion in a straight line and specular reflections 
from the walls of its container. Because they demonstrate a broad variety of behaviours (regular, chaotic [2] and 
mixed phase space dynamics) they have been readily used as models in theoretical and experimental physics [3H6] . 
They are widely applicable because their dynamics corresponds to the classical (short wavelength) limit of wave 
equations for light, sound or quantum particles in a homogeneous cavity. 




FIG. I: (Color online) a) Simple mushroom, b) Elliptic mushroom with triangular stem, c) 'Honey mushroom' with 3 integrable 
islands and 1 ergodic component. 

The mushroom billiard is constructed by a convex semi-elliptical (including semi-circular) 'hat' attached to a 
'stem' such that their intersection is smaller than the diameter of the hat. Examples are shown in Figure [T] It 
is special in that under certain conditions [3 15] , it forms a class of dynamical systems with sharply divided phase 
space which is easy to visualize and analyse. For example, the phase space of the mushroom shown in Figure [l]a) 
is composed of a single completely regular (integrable) invariant component and a single connected (topologically 
transitive) chaotic and ergodic component, in contrast with other generic mixed systems such as the standard 
map [5], where KAM hierarchical islands form a dense family in the neighborhood of each other. Interestingly, 
mushrooms can also be designed to have an arbitrary number of integrable and ergodic components (see Figure [I] 
c)). Hence, mushroom billiards are paradigmatic models for studying the phase space dynamics near the boundary 
of integrable islands. However, one must note that small perturbations (imperfections) to their boundary may 
cause the emergence of KAM islands or even complete chaos [5J [TU] . 

Mushrooms have become increasingly interesting to the quantum chaos community because of their unusually 
simple, divided classical phase space. This has facilitated for the numerical verification of Percival's conjecture 
which states that in the semiclassical limit, eigenmodes localize to one or another invariant region of phase space 
(regular or chaotic) , with occurrence in proportion to the respective phase space volumes |I11 II2j ; recently this has 
been applied to generalise the boundary term in Weyl's law }13) . Similarly, the mechanism of dynamical tunneling 
between classically isolated phase space regions has also been investigated in the context of mushrooms [T2J [T5] and 
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has been observed in microwave mushrooms [IB] . However, although their classical phase space is sharply divided, 
generic 'simple' mushrooms (Figure [l] a)) have been shown to have an infinite number of different families of 
marginally unstable periodic orbits (MUPOs), 'embedded' in the ergodic component of their phase space [TTIfTB] . 
These MUPOs live in the mushroom's hat and resemble periodic orbits of the circle. The flow close to these 
orbits, is reminiscent of that close to KAM islands [H] and thus causes the system to display the phenomenon 
of 'stickiness', where chaotic orbits stick close to regions of stability for long periods of time [20], the quantum 
analogue of which is not well understood. 

It is worth mentioning that non-sticky mushrooms have been previously constructed using elliptical hats and 
non-rectangular stems [S] . This is because each focus of the ellipse provides a sharp boundary between rotational 
and librational orbits, and may be used as the end point of the entrance to the foot. However, in such a case, some 
care is needed with the stem's length and it's base width, to ensure sufficient defocusing. In addition to this, the 
size of the opening of the stem must also ensure a bounded number of maximum possible collisions in the hat. 

In this paper we focus on classes of mushrooms with circular hats, in which stickiness is due to MUPOs. We 
express the problem of finding them as the well known problem of finding optimal rational approximations of 
a number (section II A.). This remarkable connection made with number theory allows us to introduce and 
characterize a zero measure set of control parameter values, using continued fractions, for which all MUPOs are 
completely removed (section II B.). This set, not previously discussed in the literature, corresponds to mushrooms 
with a less sticky hat, the implications of which are yet to be studied classically or quantum mechanically and 
are hoped to be useful in various applications such as directional emission in dielectric micro-cavities [211 122j . We 
obtain upper bounds for MUPO-free and finitely sticky irrational mushrooms and also give an explicit example of 
a MUPO-free mushroom billiard (section II O). 

As discussed in Ref [53], placing a small hole on a billiard's boundary allows one to 'peep' into the system's 
dynamics. The smaller the hole, the smaller the observational effect (whether quantum or classical) on the 
dynamics. We therefore 'open' the billiard with a hole and look at the survival probability P(t), given a uniform 
initial distribution of particles. By considering linear perturbations of MUPOs, we obtain exact expressions for 
the asymptotic algebraic decay of P(t). This is done for two separate cases, firstly for MUPOs in the semi- 
circular hat of the mushroom (section III A.) and then for bouncing ball orbits in the case of a rectangular stem 
(section III B.). The explicit form of these expressions in turn allows us not only to predict but also to calibrate 
the survival probability function by changing the various geometric parameters involved. Finally, the results are 
confirmed numerically (section III C.) and then discussed briefly, including higher dimensional mushrooms and 
other implications of this work (section IV). 

II. STICKINESS IN CLOSED MUSHROOMS 
A. MUSHROOMS WITH MUPOS 




FIG. 2: (Color online) Left: Any orbit intersecting the dashed red semicircle of radius r is unstable while any orbit not intersecting 
it is stable. MUPOs are the periodic orbits which intersect the semicircle while not entering the mushroom's stem. The MUPOs 
( s »j) = (4,1) and (s,j) = (5,1) are shown. Right: Phase space plot of the ergodic component for a simple mushroom using 
Birkhoff coordinates (z,sin<9), with = 0.5 and L = 0.5 computed from a random chaotic orbit followed for 10 collisions. Here, 
z S [0,7r_R + 2(R + L)) is the arc length parametrization along the billiard's boundary, increasing from zero at from the right-most 
point of the mushroom, in an anticlockwise fashion and 8 £ (— -|, ^) is the angle of incidence at each collision. 

Stickiness [2U] is due to points or regions of phase space having vanishingly small local Lvapunov[2"l] exponents, 
typically found in the vicinity of KAM elliptic islands or in the very close neighborhood of marginally unstable 
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periodic orbits (MUPOs). Sticky orbits exhibit long periods of quasi-regular behaviour and are thus also associated 
with the phenomenon of intermittency |25) . A well known and well studied family of MUPOs are the bouncing 
ball orbits present in the chaotic Bunimovich stadium billiard [26 . They form a zero measure family of period 
two orbits [27J , trapped forever between the stadium's parallel walls. Even though MUPOs do not affect the 
overall ergodicity of the system, they govern long time statistical properties of the system, such as the Poincare 
recurrence times distribution Q(t) ~ i~ 2 [28], the rate of mixing (the rate of the decay of correlations) C(t) ~ t^ 1 
[2"5H3"T] and the long time survival probability P(t) ~ £ _1 [32] given a carefully positioned 'hole'. Furthermore, the 
exponents of these power-laws appear to be a universal fingerprint of nonuniform hyperbolicity and stickiness, at 
least for one and two dimensional Hamiltonian systems with sharply divided phase space [H]. The MUPOs in the 
mushroom's hat and in the annular billiard were extensively studied by Altmann in his PhD thesis [33] and more 
briefly in |17H19j and are common in many billiards with circular arcs, the significance of which has only recently 
been realised in the context of directional emission in dielectric micro-cavities [3TJ [22] . 

MUPOs in the mushroom are best understood when introduced geometrically. The dashed red semicircle of 
radius r in Figure [2^,) corresponds to the border between the ergodic and regular component of the mushroom's 
phase space (see Figure [2]d). Any orbit intersecting this semicircle is unstable and lies in the ergodic component 
in phase space while any orbit not intersecting it exhibits regular motion and remains forever in the mushroom's 
hat. However, as shown in the figure, one can find periodic orbits in the mushroom's hat which do intersect it and 
therefore are unstable. A compact way of describing them is given by: 

in r cos , 
a »j = cos — < o < — r = Ps.j, (1) 



s R cos 



where 




if sis even, 1, if s is even, 

V - J if s is odd, [2, if s is odd. v ; 

In Eq(l), r and R are as defined in Figure U The coprime integers s and j describe periodic orbits of the circle 
billiard with angles of incidence 8 s j = \— ^ . More specifically, s is the period and j the rotation number of the 
orbit. Ra s ,j is the shortest distance from the periodic orbit (s,j) to the origin. Rfi s ,j is half the longest straight 
line passing through the origin which intersects the unfolded (along the hat's base) periodic orbit (s,j) at equal 
distances on either side. Hence, (1) guarantees that (s, j) is a MUPO and can be oriented in such a way as not 
to enter the stem while still intersecting the dashed semicircle. Let S p denote the set of periodic orbits which are 
marginally unstable for a given p = . 

It can easily be seen that a small perturbation -q with respect to the incidence angle 9 s j of a MUPO will cause 
the orbit to precess in the opposite direction and eventually enter the stem and feel the chaotic effect of the 
defocusing mechanism [34) . However, since the precessing angular velocity is proportional to the perturbation 
strength 77 which may be arbitrarily small, the orbit will behave in a quasi-periodic fashion and entry into the 
stem may take an unbounded amount of time. Hence the term 'stickiness', meaning that orbits in the immediate 
vicinity of MUPOs stick close to the regular component of phase space for long periods of time. Note however 
that although these periodic orbits are dynamically marginally unstable, they are not structurally robust against 
parameter perturbations of p. 

It is clear from (1) that the intervals (ot a j, ft s j) are shrinking quadratically with increasing s. We see this by 
rearranging (1) into: 

j I 1 / cos & \ 

> — arccos p > — arccos f- , (3) 



S 7T 7T V cos 



expanding for large s 



j j /TTCOt^X 1 / 1 \ 

where we have set d* — ^ arccos p and rearranging once more to get 

(5) 

where we have neglected the positive terms of order ~ s -4 , thus possibly losing some of the MUPOs; we give 
explicit bounds on this term in the next section. In this way the problem of finding the elements of S p is expressed 
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as the well known number theoretic problem of finding rational approximations - of $* G (0, |). However in this 
case we have a couple of complications: the approximations are one-sided and the tolerance depends both on the 
numerical value of | and the parity of s through A in (2). 

This interesting connection made here allows one to apply well known results from number theory to the present 
dynamical system and infer useful dynamical properties about it. Altmann |33| proved that for almost all p G (0, 1) 
there exist infinitely many MUPOs and hence the corresponding mushroom exhibits stickiness and all inherent 
dynamical properties mentioned above. In the following section we remove the parity dependence in the context of 
a more general mushroom model. This in turn allows us to use properties of continued fractions more effectively to 
derive a sufficient condition so that (1) has no solutions and hence destroy all MUPOs in the hat of the mushroom. 



B. MUSHROOMS WITHOUT MUPOS 



1. Generalized Mushroom 




We have seen that MUPOs in the mushroom billiard, which give it its sticky dynamical character, are directly 
related to number theory through (5). In this section we propose a generalization of the mushroom billiard studied 
in section II A., which will allow us to efficiently use properties of continued fractions without having to worry 
about the parity of MUPOs. The main result here will be to prove the existence of a zero measure set of p = 
values for which the mushroom's hat is MUPO-free. Furthermore, we shall obtain a sufficient condition which 
explicitly describes a subset of this set. 

Consider the 'elementary cell' obtained by slicing the mushroom along its vertical axis of symmetry. Then the 
period of the corresponding (s,j) orbit is sA/2. Similarly, since we are currently only interested in collisions with 
the curved segment of the billiard, we introduce the parameter a G (0, 1) which allows the mushroom to have 
circular hats of variable size. This billiard, shown in Figure [3j is in a class of billiards considered in [35], and was 
shown to have a sharply divided phase space in [5] as long as L > 0. The boundary between the two components 
is given by the dashed arc of radius r G (0, R). Notice that because the stem is triangular, there are no bouncing 
ball orbits present. 

Periodic orbits in the hat of the proposed mushroom will now have incidence angles with the curved boundary 
given by \6 q , p \ — § — ^p- for some coprime p and q, and equation (1) becomes 

rn a 93P2L 

apn r cos 

C ° S ^^i? < ^f' (6) 
Notice that there is no longer a parity dependent A. Similarly (5) becomes 

p $* got cot R^{q) 
q a 2q 2 q 2 

where •&* — — arccosp, and R2(q) is the remainder term obtained from the Taylor expansion for large q. In the 
following we bound the argument of the cotangent by ni}* and bound R2, so that for q > Q we have 

K(q,Q) 

1 



0<^--< ^ " +^, (7) 



p an cot -f- R 2 [ q ) ^ f a7Tp 



o<---< — ^H^ + ^< 

q a 2q* q z 



,2 ' 



(8) 
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where ^2(9) is bounded by 

9 / Q!7r\ 4\ p p 3 
tan 2 ' > 1 > p p 



i? 2 (g,Q) = 



(if) 9 s - 



1 



Q / 3y ^1^2 2(1-/92)1 



a 2 ir 2 \ p 3 



2 cosf 



<^) 2 (1-P 2 (1 + 

li 

greater than max (an, -^j\Jj^r^j > following from Appendix A. 



cos( 9 fy ) 2 q' 2 ' 



(9) 



COs(2f) 2 Q 2 

' it ^;us^ — ; ~ 1 J. — ij- \ x t 7— 

Q 

as obtained in Appendix A. Here, Q is a fixed number up to which (8) must be checked numerically. It must be 



2. MUPO-free Condition 

We now turn to some number theory and introduce some basic concepts. It is well known that the best rational 
approximations of a real number £ are obtained through its continued fraction representation 36 

£ = a H — j — = [ao;ax,a 2 , ...], (10) 

ai H = — 

1 a 2 + ... 

where the quantities 0,0,0,1,0,2, ■ ■ • are called 'partial quotients' and are usually taken to be positive integers. 
Irrational numbers have an infinite continued fraction representation while rationals have finite. The n th truncation 
of a continued fraction representation gives the n th 'convergent' ^ of £. Hence irrational numbers have an infinite 
number of convergents while rationals finite. Convergents are 'best approximations' to £, meaning that there is no 
other fraction with denominator smaller than B n which approximates £ better. 
For the mushroom, we would like to find values of — = [ao] ai, 02, . . .] for which 

0<P _r < m i Q lj 

q a q z 

has no solutions since this would also imply no solutions to (6). Solutions to (11), if any, are only given by the 
convergents of ^- if < K(Q,Q) < \ [37 . Obviously, if ^— is rational then there is only a finite number of 
solutions to (11). However if ^- is irrational the answer is not so simple. We focus on the convergents | = ^ of 
^ and express it in terms of them such that 

l9* Cn+lAn + An-x 



a (n+lB n + B n _i ' 

where <^ n = [a n ; a n +i, a n +2, ■ ■ ■] is the n th 'complete quotient' of — . Hence 

A n d* A n B n _\ — A n _iB n (-1)"- 1 



B„ 



Cn+l + %f ) Bl (C 



n+1 



B 2 



(12) 



(13) 



It is easy to see that if n is even, then -g 2 - — ^- < 0. Therefore, equation (11) will not have any solutions if 

K(Q,Q)<- (14) 

Cn+l + ~B^~ 

for all odd n. Since a n +i < Cn+l < o-n+i + 1 and B g~ 1 < 1, it follows that 

K{Q,Q)<^—, (15) 
p + 2 

where p = max(a2 n ), is a sufficient condition for (11) and therefore (6) not to have any solutions. The condition 
is never satisfied if a^n is unbounded. 

The set of numbers with bounded even partial quotients as derived above has zero measure [37] and has Hausdorff 
dimension one as p is unbounded as p — > 38 . As shown by Altmann |33j , a generic mushroom will have infinitely 
many MUPOs and so it will be infinitely sticky. However, we have shown here that there are infinitely many values 

of — , and therefore p, for which MUPOs in the hat are finite or completely removed. Since the smallest possible 

_ 1 

value of p is one, for the original mushroom with a = \, if p < ((^p) 2 + lj ~ 0.390683, (15) gives a sufficient 
condition for (11) not to have any solutions and therefore describes a mushroom with no MUPOs in its hat. 
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C. MUSHROOMS WITH LARGE STEMS 
1. MUPO-free Example 

For larger values, \ < K(Q,Q) < 1, solutions to (11), if any, are given by the convergents ^ and also by the 
so called 'intermediate' convergents of the form [37], where c is an integer such that 1 < c < a n+2 . 

There are however values of *&* with K (Q 7 Q) > \ satisfying (14) such that the corresponding mushrooms will 
have no MUPOs. An example of such a mushroom is p — cos ( 5 \^ nj ~ 0.64013 which has K(q, 95) < 0.6549 
and 2d* — [0; 1, 1, 3, {1, 4}] (where we have numerically checked the absence of MUPOs up to q — 95). Here, 
the odd convergents of 2d* satisfy < ^ — ^- = where K n = ^Cn+i + f° r °dd n > 3, where 

Cn+i = [1; 4, {1, 4}] = 4(1 + y/2). It is an easy exercise to show that for all odd n > 3 



B n = - ^a_A|_ — a + X: 
B„-i = ijV-Af +/3+A 5 ), 



(16) 



where A± = 3 ± 2^2, a± = 12 ± 7^/2 and /3± = ±26 + 19\/2 are all positive numbers. Hence £ 1 



2 



is strictly decreasing with n and therefore if„ is bounded by 



#(g,95) <if 5 <X n < -L, (17) 



for all odd n > 5, where ps 0.706. Similarly for the intermediate convergents of 2d* we have that 

cA n+ i + A n _ ^ _ cA n+1 + An _ C n+2 A n+ i + A n 
cB n+ i + B n cB n+ i + B n ( n+2 B n+1 + B n ' 

which for odd n > 5 simplifies to 

2 + 2V2-C _ K n (c) 



{cB n+1 + B n ) (B n+1 (2 + 2V2) + B n ) (cB n+1 + B n ) 2 ' 
since Cn+2 = [4;1,{4, 1}] =2 + 2\[2. Hence, using (16) and a similar argument as above K n (c) — 4+4 ^ c 

8(5—/2) 



(19) 



i2 - : "" i; - 2 - 2 ^ 2 ) 7 is bounded by 



K(q,95) < K 5 (l) < K n (c) < 4 + ^ C ' , (20) 

for c = 1,2,3 and odd n > 5 where -^5(1) ~ 1.237. Therefore, p = cos (^^23^ 7r ) describes a mushroom with no 
MUPOs in its chaotic region. 

2. Supremum of MUPO-free Values 



From the example above we can now use similar arguments to establish that MUPO-free values of p exist up 
to In other words sup(p G (0, 1) : S p = 0) = 4|. To see this, let K(Q,Q) denote the value of K(Q,Q) at 

p = ^j. Then from equations (8) and (9) K(Q,Q) < K(Q, Q) = § + j^s + 0(Q~ 4 ) for < p < ±. Now 
consider for m G Z + large 



^--miU>n»='^^ = l + ± + o(±), ( 21 , 



so that p = cos-7T?9* = ^= — 8v ^ m + 0(m 2 ). We first look at the odd convergents of 2d* as in (13) 



0<£-2W S §, (22) 
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where K n — 



(C«H 



Bra- 1 



we obtain that for odd n > 3 



and Cn+i = [l;m, {l,m}] = \ - 
B n -i (2-A_) + (A+-2) 



where A± 
to 



i(2 + m±v / 4m^ 
= 1-1: + A 



ro(4+m) 



(l + 2m-A_) -(l + 2x- 
3) . Thus if„ > if j 



A+) 



Via a similar manipulation as in (16) 



(23) 



l 

4 m 



3m+^/ m(4+m) 

C(m -3 ) with n and therefore if„ > K(q,B n ) 



-0(m 2 ) converges exponentially 
ps J for large enough to and n. 



Similarly, when looking at the intermediate convergents of 2d* as in (18) and (19) such that 



K n {c) 



(Cn 



+2 



(Cn+2 



where c = 1, 2, ... (to - 1), Cn+2 = [m; 1, {to, 1}] = ^ - 
that if n (l) < K n (c) and if„+i(c) < AT„(c). Hence, since if„(l) converges exponentially to 



(24) 



and n 6 " can be obtained from (23), we find 



l-2m 



m 2 — my m(4+m) — 1 



+ C(to~ 3 ) with n, then K n (c) > K n (l) > 



1— 2m 



1 > 1 — x y ' ' ' 1 w j i i .i ) / . i 1 1\ i > t v \ \ j .i \. n \ .i. i . t / ^ 

m^—mW m(4+m) — 1 

and n, thus verifying our claim above for the supremum of MUPO-free mushrooms 



> K(q,cB n+ i + B n ) S3 j for large enough to 



3. Supremum of Finitely Sticky Irrational Values 

Another claim we can make is that sup (/? : #6> p < oo, i arccosp ^ Q) = ^g 4 ^ ~ 0.7864. To see this we take 
the leading order term of K(Q 1 Q) as Q — > oo and equate it to one, so that = 1. Now since p = cos7ri9*, 

then 2d* = f arccos yj=s = [0; 2, 2, 1, 3, 1, 1, 1, . . .] < [0; 2, 2, 1, 3, 1, 2, 1, {1, to}] = 2$*. It follows that the 

corresponding Jf (Q, Q) < 1 in the limit Q — > oo. Now since we may augment the tail of the continued fraction 
expansion of 2d* as done above by the transformations a v — > a u + 1 and C^+2 [1; m, {1, m}] for any even i/, 

then K n = B% — 2d*^j will converge exponentially to some function f(m) = 1 — ^ + Ri(m) with n for some 

constant and Ri(m) = 4-^ for some < £ < to. Therefore, as f — > oo, (#* — i9*) — » + and K(Q,Q) — > 1. 
However we may always choose m and n big enough such that if n > K(B n , B n ). A similar statement can be made 
for the intermediate convergents of 2d* . For values of p > ^ 2 , K(Q,Q) > 1 and therefore all convergents of 
2d* £" Q are solutions of (11) [33] hence describing mushrooms with infinitely many MUPOs. 

In the next section we investigate the stickiness in the hat of the mushroom in the context of escape through a 
small hole placed on the stem of the mushroom. 



III. ESCAPE FROM THE MUSHROOM 



In the previous sections we have investigated the dynamics of the mushroom billiard and specifically focused at 
the chaotic region of phase space, close to the regular island. We have seen how MUPOs come into existence, how 
they affect the dynamics of orbits in their immediate vicinity but also how they can be removed. In this section 
we shall address the problem of escape through a small hole placed on the stem of the mushroom. We shall derive 
exact expressions, to leading order, for the survival probability for two specific cases of the mushroom. 

The uniform (Liouville) distribution projected onto the billiard boundary has the form {2\dQ\)~ l d.z dsinf?, where 
\dQ\ is the perimeter of the billiard while z £ (0, \dQ\] is the length parametrization round the billiard boundary 
and 9 £ (— §, §) is the angle of incidence with it. This is the most natural choice for an initial distribution of 
particles. Given such a distribution, the probability P(t) that a particle survives {i.e. does not escape through k 

small holes) in a fully chaotic billiard up to time t decays exponentially ~ e~ 7 *, with 7 ps fefragr [23! to leading 
order, where (r) is the mean free path, Ci is the length of each hole and \Q\ the area of the billiard. \dQ\ is easy to 
calculate and hence for the remainder of the paper we leave it in general form to allow for different shaped stems. 
For billiards with mixed phase space such as the mushroom, P(t) for long enough times is expected to decay as 

P(t)^A + B(e-^ + ^), (25) 
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where we have neglected terms of order t 2 . A is the measure of the integrable island given by: 

A = 4(2\dQ\)- 1 Ry/l - p 2 - pR arccos p + - p) , (26) 

and B is its complement (B = 1 — A). In (25) we have assumed that the holes are placed well in the ergodic 
component of phase space and therefore 

(27) 

while the mean free path in the ergodic component is now 



P 2 



r *(\Qs\ +R 2 arcsin p + pR 2 Jl 
c M B\dQ\ 



where \Q S \ is the area of the mushroom's stem while c„ and are the invariant probability measures of the ergodic 
component for the billiard flow and map respectively. 

Algebraic decays, of the form j in (25), originate from the stickiness exhibited, which is itself due to the nonuni- 
form hyperbolicity of intermittent systems as discussed in the previous sections. It is a geometrical description of 
the constant C that we seek here. In the case of the stadium billiard for example, near-bouncing ball orbits were 
studied and such a constant was successfully calculated in Ref [32]. In the following two subsections we attempt 
to do the same, first for the MUPOs living the mushroom's hat and then for near-bouncing ball orbits present in 
mushrooms with rectangular stems. 



A. STICKY HAT 




FIG. 4: (Color online) Left: Mushroom with triangular stem. Right: Image reconstruction trick at the base of the mushrooms hat. 
Orbits entering the lower semicircle through the thick blue horizontal line of length 2r are assumed to escape through the hole e soon 
thereafter. 



We consider a mushroom with a central triangular stem and circular hat as shown in Figure |4j hence removing 
any bouncing ball orbits between parallel walls. The only source of stickiness is hence found in the immediate 
vicinity of the MUPOs in the hat and therefore the algebraic decay j is equal to the measure (relative volume 
occupied in phase space) of the set of quasi-periodic initial conditions which do not enter the stem until a time t. 
This approximation is well justified since any orbit entering the stem has a very small probability of re-entering a 
'sticky' mode before suffering escape according to the exponential decay of (25). 

We use the image reconstruction trick [7] and neglect collisions with the base of the mushroom's hat. Hence 
the dynamics in the hat remains unchanged while a horizontal slit of length 2r centered at the origin corresponds 
to the stem's opening. We parametrize the now circular boundary by the angle (f>, where (f> £ (0, 2w) increases 
anticlockwise as shown in the right panel of Figure [4j Now, it is easy to see that each initial condition ((f), 9 s j ) is 
a MUPO if the collision coordinate (f> satisfies: 



Xs-l 
k=0 



(29) 
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with 



7T 

A 

7T 

A 



arccos (p sin( 



arccos (p shit 



. d) + (* " 1) 



2tt 

a? 



(30) 
(31) 



where p = ^ and the angles 



are taken modulo 2n. Each MUPO then defines a dashed, horizontal line in the 
b — plane (the phase space), and each dashed line has length fa — fa = ^jr — 2 arccos sin # S) j) . Notice that 
f>i and 02 are not defined if sin 9 > p. 



c) e 

0.952 



6) 
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FIG. 5: (Color online) Phase space plots of the Circle billiard with a vertical slit hole in the center (see also Figure |4jl. The dynamics 
are same as for the hat of the mushroom billiard in the sense that the ICs shown in blue will not escape through the slit and hence not 
enter the stem of the mushroom for at least N collisions with the boundary. In all plots above we have used p = 0.815 and N = 200. 
Plots b) c) and d) are magnifications of a), showing in more detail the MUPOs (4,1), (5,1), (66,13) and their surrounding sticky orbits. 
It is clear that MUPOs accumulate closer to the boundary of the integrable island at 8 = arcsin p. The white dashed lines in b) and 
c) are the analytic prediction given by (29-32).) 



To help visualize how the long surviving initial conditions near the above described MUPOs populate the phase 
space, we turn to some computer simulations. Initial conditions near the integrable island's boundary are chosen 
randomly so that <fi g (0,2w) and 9 € (0, arcsin g). The ones that survive for at least N collisions with the 
boundary are shown in the top left panel of Figure [5] for parameters N = 200 and p = 0.815. We notice that 
for the selected value of p, the most dominant MUPO is the square with (s,j) — (4,1) (see also Figure [SJa)). 
In Figure ^) one can identify the pentagon orbit (s,j) — (5, 1) which like all odd s-orbits has twice its period 
(As = 10) of surviving intervals along the horizontal line 51 . Further magnification into the phase space reveals 
the (s,j) = (66,13) orbit (see also Figure [5jl)) and then an accumulation of higher order orbits closer to the 
island's boundary at arcsin p. The next MUPO is (s,j) — (920, 181). 

We introduce a small perturbation 77 <C 1 in the angle 9 s j of each MUPO and expand (30-31) to leading order: 

Hi^., + v./.) = '.'ii>.-.,-./.)+ ! 1 - , C0S ^ J = 1 // our). W) 

p 2 — sin 2 9 s j 

/,./.•) = ..,. f j i - , oos0 °' j = 1 // rour). c-wi 

p 2 — sin 2 9 s j 
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We also impose a time constraint such that the perturbed MUPO will survive up to time t by requiring that 

cj)><t> 1 (9s,j + V,k)+2 V N, (34) 
4><Me. d + ri,k)+2TiN. (35) 

where N — 2 R C oa{e — 7 +rjj * s ^ ne num ber of collisions in time t. Expanding (34-35) to leading order defines in total 4 
lines which form a quadrilateral in phase space with area A s j which can be integrated with respect to the invariant 
measure (2\dQ\B)~ 1 d(j) dsin# to give: 

_ 8flcos 2 sd (tt - sAarccos (psm6 sJ )) 2 ( 1 \ 
Asj " 2s^\dQ\Bt W' 

to leading order in t. There are 2Xs such quadrilaterals due to ^-symmetry, however only half of the total area for 
each MUPO lies in <fi € (0, tt), which corresponds to the actual mushroom's hat. As for the initial conditions on 
the straight segments of the hat, since the billiard map is measure preserving, only 2Xj quadrilaterals are mapped 
onto them. Hence overall we obtain: 

\= £ \(s + 2j)(A sd -6 sd )+o(±) ) (37) 
where S p was defined in section II A. and 

5 sj = h^ 2 > ifco ^="> (38) 
I 0, otherwise. 

accounts for the possibility that a MUPO is situated exactly on the border of the chaotic region and therefore can 
only be perturbed from one side. The sum in (37) converges since the elements of S p , if any, are distributed with 
a bounded density with respect to In s. Also, notice that C does not depend on the size or position of the hole on 
the stem. Numerical simulations of the survival probability function are performed and discussed in section III C. 

B. STICKY STEM 

In the previous section we derived an expression to leading order for the asymptotic behaviour of P(t) (see 
equations (25-29) and (36-38)) for a sticky mushroom with a triangular stem. Here we investigate the stickiness 
introduced by the bouncing ball orbits present in mushrooms with rectangular stems of length L and a hole of 
size e on one of the two parallel segments as shown in Figure [6] A method for calculating the contribution of 
these orbits to P{t) was devised and explained in detail in [32]. Here, we follow this method and obtain an exact 
expression to leading order for the survival probability of the mushroom billiard. In doing so we discover an 
interesting discontinuous dependence of P(t) on p = ~ and also show that in the limit p — > 1 the expression for 
P(t) reduces to the one obtained in [32] for the stadium billiard. 

We first split the billiard's boundary dQ into four, non-overlapping, connected segments: dQ b s ,dQf 7 dQ\ and 
dQ c hl referring to the stem's base, the stem's parallel walls, the hat's base and the hat's curved segment respectively. 
We parametrize the right parallel wall of dQ™ by s £ (0, L) such that the interval (h~,h + ) defines the hole of 
size e as shown in Figure [6] It is now easy to see that initial conditions (ICs) (x, 9) with < x < h~ and 
< 8 < arctan cannot jump over the hole and therefore do not interact with the mushroom's hat. Such orbits 
behave in a completely regular manner and therefore can be integrated directly to give 



n / .arcsin(h /t) ph / .arcsin(2/i ft) p2h ir, i — \ 2 i 1 1, - \2 

where we have neglected terms of order ~ t~ 2 and multiplied by 2 due to the horizontal symmetry of the billiard. 
Similarly, ICs with h + < x < L and > 6 > — arctan give 

>((i-h+)/t) f L-h+ (r-h+\ 2 

1 cos6>dsd6> = 1 ; . (40) 



2\dQ\Bj J tsine 2\dQ\Bt 

ICs from dQ h s have contributions of order ~ t~ 2 to P{t) and therefore are ignored. 

As expected, the survival probability at long times is proportional to the square of the available length on either 
side of the hole. For the remainder of this section we consider ICs (a;,, 0j) such that h + < xi < L and < $i <C 1, 
and investigate how they contribute to P(t). We let n denote the number of collisions a particle experiences from 
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■f 2rntan9i) > 
We can see that 



a 



straight to straight segment before entering the hat of the mushroom, and define d\ = L — (x 
as the distance from the edge of the straight to the point of the last straight wall collision. 
n = 2rtan8- j ' wnere L'J N are the floor and ceiling functions respectively. Note that < d\ < 2rtan^. Once 
particle enters the hat of the mushroom it is advantageous to switch to coordinates suitable for the circle billiard 
map given by (4>, ip) — > {<p + -k — 2ip, ip) such that <f> is the angular collision coordinate and increases from zero 
in an anticlockwise fashion as shown in Figure |6j while ip € (—§!§) is the angle of reflection. Note that <p is 
different from what was used in section III A. Also, ip is used instead of 8 here to distinguish between collisions on 
the curved segment of the billiard boundary (dQ^) and collisions elsewhere. Once in the hat, we neglect collisions 
with the vertical base dQ b h , by using the image reconstruction trick as before. We find that the particle entering 
the hat will first collide with dQ c h at 



di 
R 



+ (l + p)di> 0, 



and its angle will be 



1> = - 



di 
R 



(41) 



(42) 



Let 0/ he the final angle obtained when the orbit re-enters the stem of the mushroom after experiencing a reflection 
process (a series of k £ Z + collisions with <9Q^) in the hat. We thus find that 



2kd x 



-{2kp + l)6i. 



(43) 



By carefully investigating the reflection process we find that k is actually restricted to only three possible scenarios 
such that k can either be equal to 1, [~] or [— 1 + 1, depending on the ICs di), which agree with the so called 
'magic numbers' from Ref [TU]. This can be seen if one looks at the least number of iterations of the circle billiard 
map before the orbit described by (41) and (42) intersects the horizontal slit hole: 



inf jeZ^ 



(2j- 1)^ + 



< P 



(44) 



In equation (45) below we have substituted the possible values of k into (43) and also calculated the values of 
0i for which each collision scenario corresponds to: 



( Ml 

R 



(2p + 1) 9i < 0, 



7r -(2Cp+l)^>0, 
i ^^-(2(C + l)p+l)^>0, 



2p(p+l)R 
di 

2pR 

(2Cp-l)di 
2(p 2 R 



< 



<6i< 



< 0,. < 



(2Cp-l)rfi 
2Cp 2 fl ' 

(2p+l)di 
2p(p+l)i?.' 



k = 1 collision 

k = C collisions 

k = (C + 1) collisions 



(45) 
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where we have set ( = \^~\ . Note that d\ is a function of both Xi and Oi. The first inequality on Oi (k — 1 collision) 
seems to suggest that Oi is unbounded, however this is not the case. This can be seen in an example situation 
plotted in Figure |7j where we have made the substitution ui 



■inOj 



£ (0,p). Notice that if p 1 is an integer, 




1— 1 




_L 

I 








II 

; 


II 

/ 



FIG. 7: (Color online) Reflection process in the mushroom's hat described in (43) using the substitution ui = 2 ^ 1 9 £ (0, p) with 
0.01. The Red, Blue and light Blue lines correspond to k = 1, f + 1 and f reflection process respectively as 



cos ^ 5+ ^?^' r and di 



p — v wo 

inscribed in the figure 



then the ^-collision processes in (45) is no longer attainable and we only have two possible collision scenarios. It 
is interesting to note that if p = 1, equation (45) reduces to equations (10-11) of [32] which refer to the stadium 
billiard's reflection process with the curved segment. 
We now formulate the time of escape for ICs (xj, Oi): 

t(x l ,9 l ,k)^ L —^ + L J - t ^+2R(p + k + l), (46) 
Vi \0f\ 

where we have taken small angle approximations, and substitute the values of Of and k for each collision scenario 
to get three equations for the time to escape. Each one of these equations describes conic sections since they are 
quadratic in both Xi and Oi variables. Rearranging to make Oi the subject, we obtain three hyperbolae in the 
Xi — Oi plane, describing the ICs that escape exactly at large times t. It is important to know the domain of 
each hyperbola. This can be obtained by substituting for the d\ variables into the inequalities of (45), and then 
rearranging for Oi. These inequalities are given below for the corresponding collision scenarios: 

> „■ • L-Xi (2Cp- l)(L-Xi) 

k = ( collisions , — — < Oi < y SP 



2p(l + n)R 2p(pR(l + 2n) - 2pnR ' 

k = (c + 1) collisions , J'lrTr: » < *« < {2p +1)iL ~ Xi) 



2p(pR(l + 2n) - 2pnR 1 2pR(p + 2np + 1 + n) ' 

k - 1 collision , (2p+l)(L- Xi ) 

2pR(p + 2np + I + n) V ; 

For n = 0, 1, . . . and for t large, we plot the three hyperbolae from (46) subject to (47) and the three straight 
lines from (45) onto the Xi — 0i plane (see Figure [8]). These define an area in phase space which corresponds to 
the ICs that survive at least until time t for fixed n. The various colors indicate the type of reflection process k 
the ICs experience in consistence with the ones in Figure [7] Notice that as the number of collisions n with the 
straight segments increases, the area of interest tilts and stretches in a non-overlapping fashion. To obtain the 
contribution to P(t) of these long surviving ICs, we must integrate each non-overlapping area and sum them all 
up. Note that the invariant measure will be assumed to be dp, = (2\dQ\A4)~ 1 dO i dxi here since Oi is small and 
thus dsin^w d#i. 

The corners of each enclosed area A — G, as shown in Figure [H| for each value of n are given in Appendix 

B. There are various issues which one needs to consider in order to obtain correct asymptotic expressions for 
the areas. Firstly, one needs to approximate all the hyperbolae by straight lines. This is done by joining the 
corners A — G and thus forming an irregular polygon. For example, the hyperbola between A and F, which comes 
from t(xi, Oi, 1), is approximated by a straight line joining A and F. Similarly, for the hyperbola joining B and 

C, which comes from and for the hyperbola joining D and E which comes from t(xi,0iX + 1)- The 
remaining edges are already straight lines and thus need no approximating. As shown in |32j . the error in these 
approximations is O (t~ 2 ) and hence meets our required asymptotic accuracy. 
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FIG. 8: (Color online) Area enclosed by equations (46) subject to (47) and equations (45) in the Xi — 9i phase space for n = and 1, 
using p = 0.6, L = 1 and t = 50. The colours used are in consistence with the ones in Figure The dotted, dashed and solid black 
straight lines come from the inequalities in equation (45). The corners A — G are defined in Appendix B. and are each highlighted by 
a black dot for n = 0. The dashed vertical line at Xi = h+ shows how the hole truncates the area of interest. The area defined for all 
n, corresponds to the ICs that survive at least until time f = 50. 



Another issue to be dealt with is the position of the hole which restricts the irregular polygons in Xi E L). 
This forces a deformation by truncating each polygon from the left each time one of its corners surpasses the 
hole's position as seen for example in Figure [8] This is due to the tilting effect caused as n is increased. 
Following i32] again, we expect 7 different sums since there are 6 corners (A — F), each of which will intersect 
the hole at h + at different values of n. We thus solve for n and find that the leftmost corner A T , = h + when 



n A 



E (l+2p)-2-20p-12p -4p : 



4p+i2p 2 +8p 3 J ' Similar expressions have been obtained for all other corners (B — F) and are 

given in Appendix C. Interestingly, we find that the order in which the corners A— F coincide with the hole's position 




FIG. 9: (Color online) n^ — np are defined in the text above and given in Appendix D. The figure shows how they vary discontinuously 
as a function of p = £ (0, 1) for t = 10 4 . 



depends on the system's control parameter p — Their order alternates between ua < tib < no < tie < nc < tif 
and Ha < tid < n>B < n c < n E < nF for p € (0, 1), which is shown in Figure [9j This is due to the discontinuity 
introduced by the ceiling function in £ for hb — He, hence the lower and upper bounds of the 7 sums will depend 
on the above order, and so will their arguments. Altogether we write: 

n B n D n E nc n F oo 

£A+ £ P2+ £ Ps+ £ iU £ h+ £ h+ £ A, (48) 

n— n— n.4 + 1 n— + 1 n—njj-\-l n— + 1 n— nc + 1 n— np+l 
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tia n D n B no »E np 



£A + E ft+ E ^+ E ^+ E p s+ E fl>+ E ^ ( 49 ) 

n— n— + 1 n-njj+1 n-tig+1 n-nc + 1 n— n— + 1 

where • and " are used to distinguish between the two orderings described above. Pi and Pi, i = 1, ... 7, are the 
respective areas of the polygons which we are summing over. Note that Pi = Pi and Pj — Pj. The process of 
finding all the P, and Pi is long but fairly elementary. 

We now obtain leading order expressions for each sum in t. The way to do this is similar as in [32], where a 
more detailed explanation of the method can be found. First we substitute t = h, and then n = - into the p 
and the Pj, such that u is small and v = O (1). We Taylor expand Pi and Pi into series up to order u 2 and then 
reverse the substitution by setting v — nu, thus effectively incorporating the large n into the leading order term of 
each series expansion. Now each sum can be simplified into expressions involving polygamma functions of order 
and 1. 

The polygamma function of order i is defined as the (i + l)th derivative of the logarithm of the gamma function: 



cl( 2+1 ) 

* (z)= d^) lnr(z) - (50) 



The polygamma functions are of the form z = £- +c, where a, b and c are real constants, and can thus be expanded 
as a Taylor series to leading order as follows: 

* (0) (£ +c ) =ln (l£ i l) +o(w) ' (51) 

*(*>i) (~ + c) = (-l)C*-i)(i - l)l(— )* + 0{u l+1 ). (52) 

With these approximations at hand, we obtain expressions for the sums in (48) and (49). We only present here 
the first of the approximated sums and include the rest in Appendix D: 

n A n A , + , 2 

EA = £A = !^> (53) 



where we have neglected terms of order ~ t . Altogether (48) and (49) take the form 



(L-h+f 
4C(1 + C)pt 



.2 <r„„3j_ c . A 



(2p+l)(2Cp-l) 2 



+ ln((2p+l) J1 (2Cp-l) 32 ) 



(54) 



where the coefficients Si (i = 1 . . . 4) and jj (j = 1,2) are given in Appendix E for both orderings : and ~. It 
remains to multiply (50) by 2 due to the horizontal symmetry of the mushroom, and normalize by 2|9<5|S to 
obtain a probability. The sum of expressions (39-40) and (54) depending on the value of C, therefore gives the 
asymptotic contribution of the long surviving near-bouncing ball orbits j to the mushroom's survival probability 

Interestingly yet reassuringly, in the limit of p — > 1, £ = 2 and the complicated expression for (54) reduces 
to ( L ~ h l ^ 31n3 + 2 ) ^ wn i cn j s exactly what one would expect since in this limit, the mushroom billiard is reduced 
to the half-stadium billiard 32 . In the opposite limit where p — > 0, C — ^ oo the mushroom's stem shrinks and 
expression (54) has asymptotic expansions of 



- t Q(i - h+f - 2(1 + Q(L - h+fp + 3((L - h 



+ ) 2 p 2 \+0(r 6 ) (55) 



and 

1 

7 



~(L - h+) 2 + 4(1 + 0(L - h+fp - 4(2C - 3)(i - h+j 2 p 2 ^j + 0(r 3 ), (56) 

for the two orderings • and • respectively, indicating that the discontinuous dependence on C persists, hence this 
limit is in some sense ill-defined. 



NUMERICAL SIMULATIONS 



Having obtained exact leading order analytic expressions for all the parameters appearing in (25) we now 
numerically test their validity by plotting the conditional probability P e (t) that a particle survives up to time t 



given that the particle is chosen uniformly from the ergodic component of the billiard flow (see Figure 10). 



P e (t) = (P(t) A)/B = e-V + C - + o(^\. (57) 
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FIG. 10: (Color online) Numerical simulations of P e (t) defined in (57) are plotted on a logarithmic scale using 10 8 chaotic 
ICs as a function of t. The parameters (r,R,L,e) used for the triangular stem (Left) are (cos 0.34847T, 1, 1, 0.048) such that 

S p = {(20, 7), (66, 23), (376, 131)}, while for the rectangular stem (Middle) (cos(^±3^7r) ,1,1,0.02), with h+ = 0.3 such that the 
mushroom's hat has no MUPOs. The MUPO-free mushroom (Right) has parameters (cos 0.3484-7T, 1, 1, 0.0371). The blue curves are 
the analytic predictions while the numerical data correspond to the empty circles. The insets are plots of tP e (t) showing the agreement 
with the analytic expressions for ? . 



The plots are purposely chosen (from many more) to portray and verify the results obtained in the present paper. 
Three different mushrooms are simulated: one with MUPOs present only in the hat {Left), one with bouncing ball 
orbits in the stem and a MUPO-free hat (Middle), and one with no MUPOs at all (Right). Different hole sizes give 
different exponential escape rates 7 and in turn cross-over times to a power law decay. The empty black circles in 
the plots correspond to the numerical data while the blue curves give the analytic predictions of (57). Although 
each simulation consists of 10 8 chaotic ICs, we were not able to pick up any power law decay in the MUPO-free 
mushroom (Right) and hence any clues of stickiness. 

IV. CONCLUSIONS 

In this paper, we have attempted to quantify the stickiness observed in the mushroom billiard by placing a hole 
in its ergodic component and looking at the survival probability function P(t) at long times (see eq (25)). Our 
analytic predictions are in good agreement with the numerical simulations performed and therefore confirm that 
P(t) ~ j for long enough times. Also, their good agreement with the constants C derived in sections III A. and 
B. for MUPOs present in the hat and in the stem respectively, implies that these MUPOs are indeed the primary 
causes of the power-law decay. This observation in turn applies to the Poincare recurrence times distribution Q(t) 
studied in Ref [25J and the rate of mixing of the ergodic component [3T] . 

The explicit expressions obtained here for C, allow one not only to predict but also to calibrate the asymptotic 
behaviour of P(t). Also, we have shown that these distributions as well as the overall existence of MUPOs in the 
hat are sensitive to the system's control parameter p — whilst only the near-bouncing ball orbits' contribution 
to P(t) depends on the hole's position and size. The reason for this is that the hole intersects the sticky region 
generated in phase space by the period-2 bouncing ball orbits. This creates a fictitious, time dependent 'island of 
stability' in the mushroom's ergodic component. Although orbits in it are unstable, they only experience up to one 
non-linear collision process before escaping, thus allowing us to approximate their occupancy in phase space with 
polygonal 'spikes' which we could then integrate over. In the case of the MUPOs in the mushroom's hat, we could 
easily bound the long surviving orbits by assuming that they will escape exponentially fast once in the stem. It 
is expected that the methods used here can be further generalised and applied to other mushrooms with elliptical 
hats for instance, or even to other billiards such as the annular or drive-belt stadium billiards where polygonal 
type MUPOs act as scaffolding for sticky orbits to cling onto. 

A major result of this paper is the introduction of a zero measure set which describes MUPO-free mushrooms 
(see section II B.) which to the best of our knowledge posses the simplest mixed phase space in two dimensions. 
The interesting connection between mushrooms and the diophantine approximation discovered by Altmann |33j 
cannot be directly exploited due to a sensitive parity dependence of the periodic orbits of the mushroom. We have 
overcome this complication by considering a generalized mushroom with a variable sized hat and triangular stem 
(see Figure [3]). This allowed us to efficiently use properties of continued fractions and characterize a subset of 
the infinitely many MUPO-free mushrooms. We thus obtained upper bounds for MUPO-free and finitely sticky 
irrational mushrooms and also gave an explicit example of a MUPO-free mushroom billiard (see section II C). 

Unlike the non-sticky elliptical mushrooms mentioned in the introduction, the MUPO-free mushroom with a 
circular type hat exhibits a reduced amount of stickiness (larger scaling exponent). This is due to the difficulty 
in 'finding' the foot of the mushroom by orbits which are just inside the dashed circular arc of radius r. Hence, 
the MUPO-free mushroom, although not generic, is mathematically important as it should display a polynomial 
decay of correlations of order ~ t~ 2 . 
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The results of section II also apply in the case of a three dimensional mushroom billiard with a hemispherical 
hat of radius R, a cylindrical stem of radius r and height h > and a cuboidal pedestal of base length I > 2r to 
break angular momentum conservation (see Figure 11 ). This is because orbits inside a three-dimensional spherical 
billiard always lie in the same two-dimensional plane which also contains the center of the corresponding sphere. 
The remaining (zero-measure) MUPOs in such a system are of the bouncing ball type and are found both in 
stem and pedestal. We conjecture that the corresponding mushroom has a sharply divided phase space with a 
MUPO-free hat [35]. 

Finally, one would expect to see the classical dynamical features caused by the removal of MUPOs to appear in 
the analogous quantum system in accordance with Bohr's correspondence principle. Two obvious candidates for 
such an observation are the localization of wave-functions (scars) [12] and the effect on the dynamical tunneling 
rates |14| , thus providing a new model for quantum chaos |40j as well as experimentalists with applications which 
range from semiconductor nano-structures to dielectric micro-cavities. 
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Appendix A: A Bound on R2 

To obtain a bound on R 2 in eq (8) we shall use the remainder term from Taylor's theorem several times. Taylor's 
theorem states that if / is a function which is n times differentiable on the closed interval [a, x] and n + 1 times 
differentiable on the open interval (a, x), then 

f(x) = /(a) + ^(x - a) + ^(x - a ) 2 + ... + ^-^(x - a) n + R n (x), 
1! 21 n\ 

where R n {x) — ^qr[jr ( x ~ a )" +1 f° r somc a < £ < x. 

Let c = cos (^air^j < cos7r?9*, and e — c ( — — 1^ such that eq (6) corresponds to 

c < p < c + e. 

We have to compute 
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The Taylor expansion of arccos (c + e) at e — is 



s 

arccos (c + e) = arccos(c) — ^ | + A x (AI) 
with Ai = ^±5 r£ 2 < q an( j o < a; < e. Now, for q — > oo, we have 

2(l-(c+j;) 2 )2 

ell \ a 2 ir 2 , , . 

c=(c^(f)- 1 J=^ +S2 = Bo (A2) 

with q Bo — — i m — > 2 and <7±>o = — i ~vr for some < z,u> < A . Since for q > Q we have 

f < — ^ — r < — r — v and < sin(a7rz) < sin ( ^ ) . Hence 

cos(anz) — C os(^) v ; — \ Q J 

-.Mil' — 1 " * ^ '' X i ' ' ' ' — 1 " 



Sin (^) 2 5 \ 3 3 . (an\ /sin(^) 2 5 \ a M 



and similarly 

< qB a < 



cos(^)V 



This, together with (|A1|) and (|A2|), gives 

2^1^ q 2 



Q; 2 C7f 2 X 

arccos(c + s) = arccos(c) — + C\ + Ai 



where Ci = — ?== < is bounded by 

vi— c 

. /sM^) 2 5 \ a 4 7r 4 c ^m(^) 2 5 ^ « 4 tt 4 costhT 



^cos(^) 4 6cos(^) 2 y / g 4 71^2 - \ v cos(^) 4 6cos(^) 2 y / g 4 ^1 - (costh?*) 2 
and < by 

I A | < gO+jo) 2 B 2 < f, , V 2 * 2 \ «W 

\Ai\ < r c #o S 1 H /a7r\ 2 ^, 2 a — 

2(l-c 2 (f + B ) 2 ) 5 V cos (l^) Q J 2cos(5f) 4 (l-c 2 (l + 5 ) 2 ) 5 <7 4 



" I 1+ cos(^) 2 Q 2 



a 4 7r 4 (cos7r^) 3 

3 

2c os(^)4 (i - (cos7r^) 2 (f + cos ^; v ) 2 )' g 4 



for a > — a / cos . In the same way 

a 2 C7T 2 1 C£ 2 7T 2 COS 7T$* f 



2 V / T^ 2 g 2 2^/T^ 2 g 2 
with C2 < and bounded by 

a 2 7T 2 £ f a 2 7T 2 S COS7TI?* 1 Q 4 7T 4 COS7ri? 



+ C 2 , (A3) 



Co < = ^ < , " „ < 



2yi _ ^ 2 g 2 " 2 v /f-(cos7r^*) 2 g 2 " 2cos(^) 2 v/f - (costtz?*) 2 g 4 
Finally we must bound 

f 1 f 



Vl -C 2 y/l - (COSTTT?*) 2 + ((C0S7TZ?*) 2 - C 2 ) V /f-(cOS7Tz9*) 2 + ^' 

where we wrote = (cos7r$*) 2 — c 2 so that we may expand for v small 

1 



V /l-(c0S7Tt9*) 2 +^ y/1 - (COSTH?*) 5 



where Da = T , with < y < v and so it is bounded by 

2(l-(cos7ri9*) 2 +y) 2 

.„ . v 2ecos7n9* a 2 7r 2 (cos7n?*) 2 

Dn < — r < r < V '- 



Therefore for ( A3 1 we have 



2k3 2(1- (cos th?*) 2 ) 3 cos(^) 2 (l- (cosTT^*) 2 )' 9 2 ' 



a 2 ir 2 c 1 a 2 7r 2 cos7n9* 1 

^2 + 



where 



2\/r^c 2 g 2 2^1 - (cos to?*) 2 q 
a 4 7r 4 (cos7r?9*) 3 



|C 3 | < 



2cos(^) 2 (l-(cos7n?*) 2 )5 9 4 ' 



Putting everything together, for q > max (<3, j^if%) \ZtT^^T ) we have 



a 2 c7r 2 1 



arccos(c + e) = arccos(c) + Ci + 

2vl - c 2 <r 



2 7r 2 cot ird* 
~2q~ 2 



arccos(c) ^ h Ci + A a + C 2 + C 3 , 



where the remainders have magnitudes bounded by 

/ sin(^) 2 5 \ aMcotTnT 



cos(^) 4 6cos(^) 2 J q 4 

a 2 n 2 \ a 4 7r 4 (cos7rr) 3 



\Ar\< |1+ „„^ 2Q2 a , 

® ' 2cos(^) 4 (l-(cos7r^) 2 (l+ cos g; v ) 2 ) 9 4 



cos( £ 
aV cot 7T1?* 



- 2cos(^)V 

aVfcotrff 
1 31 " 2cos(^)V ' 

and are all negative. 
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Appendix B: Corners of the polygonal Area 

The corners of the polygons (for fixed n) as shown in Figure [8] are found by solving for the intersections of the 
various curves and lines obtained from equations (45) and (46): 



Be, 

c 6i 

De, 

E Xi 
E 6i 
Fx, 

Fe, 

Gxi 



-2{h+p{\ + 2p)(l + n + p + 2np) - L{\ + p)(l + 2(2 + n)p + (2 + 4n)p 2 )) + R(-2LpR - LR)t 



2 + 2(4 + n)p + (6 + An)p 2 + R(-2pR - R)t 
(h+ -L)(l + 2p) 2 



2R(1 + p(4 + n + 3p + 2np)) - t(l + 2p) ' 

-2h+Rp(l + 2p)(l + n + p + 2np) + £(-2(1 + ()R + t + 4(1 + 3( + 2(1 + Qn)Rp 3 ) 
{2C,p - l)(-t(l + 2p) + 2R(1 + C + 2(p + p(4 + n + 3p + 2np))) 
2p{t - C,t{\ + 2Lp) + fl(-3 + 2C 2 + 2L{n + C(4 + 2C + n))p)) 
+ (2(p - l)(-t(l + 2p) + 2i?(l + C + 2Cp + p(4 + n + 3p + 2np))) ' 

(h+-L)(2p+l) 2 

(2(p - 1)(-£(1 + 2p) + 2R(l + C + 2(p + p(4 + n + 3p + 2np))) ' 
2(h+ - L)Rp(2(p - l)(Cp + n(2Cp - 1)) 



(1 + 2p){t - 2(tp + 2R(-1 - C + (-1 + 2C(1 + C) - n)p + C(3 + 2n)p 2 )) ' 

(h+-L)(l-2(p) 2 

(1 + 2p)(t - 2(tp + 2R(-1 - C + (-1 + 2C(1 + C) - n)p + C(3 + 2n)p 2 )) ' 

2{L + h+n)Rp -Lt + 4( 2 Rp((h+ - L)(l + 2n)p 2 - L) + 2((LR + Up - (h+ - 2L(n - 1) 



(Bl) 
(B2) 

(B3) 

(B4) 

(B5) 

(B6) 
4h+n)Rp 2 ) 



-t - 4( 2 Rp + 2(1 + n)Rp + 2C,(R + tp-(3 + 2n)Rp 2 



(h+ -L)(l-2(pf 



-2(R + t - 2((1 - 2C 2 + n)R + (t)p + 2((3 + 2n)Rp 2 ' 

Lt + 4( 2 LRp - 2(h+ + L + h+n)Rp + 2(L(R(-1 + 2(2 + n)p 2 ) - tp) 



(2{p-l)(2R(C + (2 + n)p)-t) 
(h+-L) 



(-l + 2(p)(-t + 2R{{+(2 + n)p)y 

L(2R{\ + p)(l + 2(1 + n)p) - t{\ + 2p)) - 2h+nRp 

(l + 2p){2R(l + p + np) -t) ' 

(h+-L) 



(l + 2p)(-t + 2R(l + p + np))' 

L, 

0, 



(B7) 
(B8) 

(B9) 

(BIO) 

(BH) 

(B12) 

(B13) 
(B14) 



where £ = ["— 1 and p 



Appendix C: Values of n when corners hit the hole 

The upper and lower limits of the sums in expressions (48-49) are the solutions for n when the Xi coordinate of 
the corners A — F exceeds h + : 



n B = 



nc 



n D 



n E 



Up 



(2 + C - |) + (6 - 4C 2 + ^(2C - 2)) p + (-16C - 8C 2 + 4C^)p 2 + (-4 - 12Qp 3 

4(l + C)(2p + l)p 2 

(2 + C - |) + (6 - 4C 2 + |(2C - 2)) p + (4 - 12C - 8( 2 + 4C^)p 2 + (-12C - 4C 2 )p 3 

4(l + C)(2Cp-l)p 2 

2C - | + (2 - 4C 2 + 2(±)p - 4Cp 2 - 4C 2 p 3 

4C(2Cp-l)p 2 
2C-| + (2-4C 2 + 2C|)p-8Cp 2 



-2+^ + (-6+^)p-4p 2 



4Cp 2 
TO 



4p(p+l) 



(CI) 

(C2) 

(C3) 
(C4) 
(C5) 



20 



where £ = [^] and p = ^. 

Appendix D: Leading order approximations of Sums in eq (48-49) 



E A = 

n— nA+l 
no 

E A 

n— n_e + l 

E A 

n— ri£> + l 

E A 

n— 

E A 

n-nc+1 



E A 

n=riF-\-l 

E A 



n— nA+1 
n B 

E * 

n— ri£> + l 

E * 
E h 

n-nc+l 

E * 

oo 

E ^ 



(fr + -£) 2 

4pt 

4(1 + C)ip 
(fc + -^) 2 

4C(i + Qtp 



(4 + 2Qp+(-8C-2C 2 )p 2 + 4C 2 p 3 
(2Cp-l)(2p+l) 



+ ln(2Cp-l) 



2p(-l + 2C(Cp) - l)(2p + C(~l + g + KP - 2p 2 )) 
(l + 2p)(l-2Cp) 2 

-2p(l + C - C 2 + 2C 3 P)(1 + P (l + 2C(-1 + Cp))) 
(l + 2p)(l-2Cp) 2 



+ (2 + C) In 



2,9+1 



(1 + 3C + C 2 )ln(2p-1) 



(V -Lf [2p(-l + 2C(-l + Cp))(2p + C(-l + P + 2Cp-2p 2 )) | | f 2 /)+ ] 



4(1 + C)tp 

(^ + -^) 2 
Apt 



(l + 2p)(l-2Cp) 2 



(2CP-1) 2 . 



2p(Cp-l)(C(2p-l)-2) 
(l + 2p)(2Cp-l) 



+ ln(2Cp-l) 



2 = 



(fr+-L) 2 (p+l) 
(2p+l)t 

- L) 2 [ (2 - 2C)p + (2 - 4C + 2C 2 )p 2 _ + ln ( 2, , - I 



Apt 

(h + -Lf 
4(tp 

(h + -L) 
4C(1 + C)ip 

(h+-Lf 



(2Cp-l)(2p+l) 



2Cp-l 



2p(l + p + 2C 2 p - <(i + pW- - 2C(Cp - i)) 

(l + 2p)(l-2Cp) 2 



2 r 



4(tp 

f -i 
4pi 



4(1 + pg(i - (C - i)p)(i + C - C 2 + 2C 3 p) 

(l + 2p)(l-2Cp) 2 

2p(l + p + 2C 2 p - C(l + p))(l - 2C(CP - 1)) 
(l + 2p)(l-2Cp) 2 



+ (1 + C)ln 



+ (1 + 3C + C 2 )ln 



+ (1 + C)ln 





1) 2 \1 


\ (2p4 


1) )\ 



(2CP^1) 2 
(2p+l) 2 

2p+l 



{h+-Lf [( 2-2Qp+(2-4C + 2C 2 )p 2 / 2p+l 
(l + 2p)(2Cp-l) \Kp~l 



(h+-Lf{p+l) 
(2p+l)t 



n— riF + 1 

where £ = ["^] and p = 



(Dl) 
, (D2) 
. (D3) 
,(D4) 
(D5) 
(D6) 
(D7) 
(D8) 

. (D9) 
(D10) 
(Dll) 
(D12) 



Appendix E: Coefficients of eq (54) 



The coefficients of equation (54) for the two different orderings : and ~ of ua — rip are 

£i = -2 - 4C + 4C 3 £i=8 + 18C + 2C 2 - 8C 3 (El) 

e 2 = -2(1 + C + C 2 + 2C 3 + 6C 4 ) £ 2 = 8 + 12C - 20C 2 + 24C 4 (E2) 

e 3 = -12C 2 - 8C 3 - 4C 4 + 8C 5 £ 3 = -16C 2 + 8C 3 + 8C 4 -16C 5 (E3) 

£ 4 = 16C 3 + 16C 4 + 8C 5 £ 4 = 16(C 3 + C 4 ) (E4) 

h = 1 + 7C + 3C 2 ji = -4 - 8C - 2C 2 (E5) 

32 = -6C - 2C 2 h = 6 + 12C + 4C 2 (E6) 
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where C = [f ] • 
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